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Abstract. The Pseudo-Nambu-Goldstone Boson (PNGB) potential, defined through the 
amplitude and width / of its characteristic potential V{(j)) = + cos{4>/f)], is one of 

the best-suited models for the study of thawing quintessence. We analyse its present obser¬ 
vational constraints by direct numerical solution of the scalar field equation of motion. Ob¬ 
servational bounds are obtained using Supernovae data, cosmic microwave background tem¬ 
perature, polarization and lensing data from Planck, direct Hubble constant constraints, and 
baryon acoustic oscillations data. We find the parameter ranges for which PNGB quintessence 
gives a viable theory for dark energy. This exact approach is contrasted with the use of an 
approximate equation-of-state parametrization for thawing theories. We also discuss other 
possible parameterization choices, as well as commenting on the accuracy of the constraints 
imposed by Planck alone. Overall our analysis highlights a significant prior dependence to 
the outcome coming from the choice of modelling methodology, which current data are not 
sufficient to override. 
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1 Introduction 

The nature of the Universe’s late-time acceleration, originally discovered via measurements of 
51 type la Supernovae [1, 2], remains unexplained. Under Einstein’s general relativity scheme, 
astronomical observations [3, 4] indicate that 68% of the energy density of the Universe 
corresponds to dark energy, a material component with negative pressure and the originator 
of the poorly-understood accelerated cosmological expansion at present. Reference [5] is a 
comprehensive review on the subject, providing a good number of candidate dark energy 
models. 

A plausible choice for explaining the characteristics of dark energy is the cosmological 
constant A, appealing for its simple equation of state w = —1 and its accurate agreement 
with observations. However, this choice has awkward implications, in particular the 120 
orders of magnitude difference between the expected theoretical value of the vacuum energy 
density and the observational magnitude of A, which raises concern. The so-called cosmic 
coincidence problem, in reference to the circumstance that both the matter and dark energy 
components of the Universe are of the same order of magnitude at the current epoch, despite 
having been so different at every other stage of the cosmos’s history, adds more doubts to the 
underpinning of the standard cosmological model. The increasing amount of available data 
and the aforementioned theoretical dilemmas have promoted the discussion of alternative 
possibilities to the cosmological constant in order to deepen our understanding of the present 
and future state and composition of the Universe. One of those possibilities is an evolving 
dark energy equation of state. 

A canonical scalar field which varies slowly along a potential V{4>), known as quint¬ 
essence, can lead to the very similar observational results as the cosmological constant, with 
the advantage of a broader phenomenology and the possibility of a link to fundamental 
physics models. Reference [6] reviews the cosmological dynamics of quintessence, providing 
approximate analytical solutions wherever possible. These theories can be approximately 
divided into two classes: tracking type, in which the equation of state decreases towards 
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It; = —1 as accelerated expansion commences, and thawing type, in which the field is initially 
frozen by Hubble friction during radiation and matter domination until late times when it is 
finally allowed to roll down the potential. 

In this paper we analyze one of the best candidates for thawing quintessence, the Pseudo 
Nambu-Goldstone Boson (PNGB) model [7]. The PNGB’s characteristic potential is suffi¬ 
ciently flat as to give an appropriately-negative equation of state with w close to —1. Previous 
efforts to constrain PNGB quintessence may be found in Refs. [8-11], while Ref. [12] carried 
out forecasts on the ability of future experiments to constrain this model. An alternative 
approach taken into the study of thawing quintessence is via a suitable approximate equation 
of state, established in Refs. [13, 14], where the parametrization is done via the present equa¬ 
tion of state Wo, the field density parameter, and the gradient of the selected potential. A 
parameter estimation analysis using this method has previously been carried out in Ref. [15] . 
The PNGB potential constraints have also been thoroughly analysed in the context of dark 
matter densities in Ref. [16]. 

This paper is organized as follows. In Sec. 2 we present the background evolution 
equations of PNGB quintessence. In Sec. 3 we describe our analysis by direct parametrization 
of the width and amplitude of the potential, performed with GosmoMG version July 2015 and 
the quintessence module included within the GAMB code [17]. We mention the prior range 
for the model parameters and the datasets utilised. Section 4 contains the observational 
bounds obtained for this model. The accuracy of the approximate approach developed in 
Refs. [13, 14] is analysed in Sec. 5 with our full data combination. The performance of Planck 
data against the Planck, JLA, HST and BAG combination is discussed in Sec. 6.1, followed 
by a discussion of the choice of priors in this particular model in Sec. 6.2 and a comparison 
against previous studies of the same potential in Sec. 6.3. We finish in Sec. 7 with some 
concluding remarks on this work. 


2 PNGB quintessence 


We assume a spatially-homogeneous quintessence field described by the scalar cj) and its 
potential The action is 


5 = y d!^Xy/^ 




+ 'S'n 


( 2 . 1 ) 


where g is the determinant of the metric g^,^, R is the Ricci scalar, Mpi is the reduced Planck 
mass, and Sm is the action for non-relativistic matter. 

For a statistically spatially homogeneous and isotropic universe with Friedmann-Robert- 
son-Walker (FRW) metric, ds'^ = —dt^ + a?‘{t)dx^, the background evolution is given by 


' -|- 3H<j) -\——— — 0 . 
d(p 

= P(f) + Pm ■ 
Pm T 3Hpm = 0 . 


( 2 . 2 ) 

(2.3) 

(2.4) 


where H = d/a is the Hubble parameter and a dot stands for a derivative with respect 
to physical time t. The pressure and energy density of the scalar field are given by P^ = 
02/2 — 1/(0) and P(j, = 0^/2 -|- 1/(0) respectively. The quintessence field equation of state is 
defined as w = Pfj,/pfi). The density parameter is = p^/pcrit, where pcrit is the critical 
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Figure 1: Equation of state versus scale factor for the PNGB potential with / = 1.4Mpi, 
for models leading to present values tco = —0.99 (red), wq = —0.93 (blue) and wq = —0.6 
(green). In thawing quintessence, w > —1. 


density of the Universe. The model under consideration is of the thawing type; which means 
it has a potential able to mimic a nearly-frozen field during the matter-dominated era, caused 
by Hubble friction, implying w near —1 at early times. Its kinetic energy contribution must 
be kept small, which means a small mass is required as well as a nearly-flat potential. The 
PNGB theory fits all of these conditions. Its potential is 


!/(</)) = 


1 + cos 



(2.5) 


where is the amplitude of the potential and / determines the width of the function. Fig¬ 
ure 1 shows some sample evolutions of the equation of state obtained numerically for models 
giving present density parameter = 0.68. The dynamics and motivational background of 
the model are analyzed in detail in Ref. [7]. 

The parameters that determine the cosmological evolution of this model are its nor¬ 
malization M^, its width /, and the initial conditions and (pi. The rapid early expansion 
of the Universe leads to sharply-decaying field velocity at early times, enabling us to assign 
(p = 0 as the initial condition for our numerical evolution. Requiring that at present the 
quintessence field has a particular density parameter allows us to use the density parame¬ 
ter as a variable and treat (pi as a derived parameter. The background evolution is calculated 
from a scale factor a = 10“®. 


3 PNGB potential analysis with CosmoMC 

The exploration of the parameter space for the PNGB potential was made using the Cos- 
moMG program, version July 2015 [19]. Within the GAMB program the quintessence module 
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Parameter 

Prior range 

f /Mp\ 

[0.1, 2] 


[0.25, 2] 


[0.005, 0.1] 


[0.001, 0.99] 

ns 

[0.8,1.2] 

T 

[0.01, 0.8] 

1OO0MC 

[0.5, 10] 

ln(10i°As) 

[2, 4] 

Wo 







Table 1: Prior ranges for cosmological and PNGB model parameters, the prior being 
uniform in the parameter quoted. The meaning of the cosmological parameters is as in the 
Planck collaboration papers [3, 18]. The final three parameters listed are derived from the 
others and inherit non-uniform priors from their relation to them. 


was updated by us to be compatible with the CAMB/CosmoMC version of July 2015. Per¬ 
turbations in the quintessence field are fully taken into account for the calculation of the 
cosmological observables. Additionally, a modification of the original code was implemented 
to solve a starting point issue: in order for the evolution to commence, the amplitnde of 
the potential has to be large enough to allow the quintessence density parameter 11^ to have 
a valne corresponding to the observed dark energy density. The random search nature of 
the MCMC code (at least at the beginning of the parameter space exploration) caused the 
program to stop through failing to meet this condition, therefore not allowing the code to 
calculate a likelihood. Instead of stopping the code after an unsuccessful initial setting, the 
unsuitable parameters were assigned an improbable negative logarithmic likelihood of 10^*^ 
(the standard valne for the program to deem a set of parameters unlikely) therefore rejecting 
them but allowing the rest of the estimations to continue. 

The choice of prior ranges for the standard cosmological parameters were taken as in 
the Planck collaboration 2015 analysis [3]. These, as well as the added free components that 
characterize the PNGB model, are displayed in Table 1. Additionally we impose a range on 
Hq of 50kms“^ Mpc“^ < Hq < 80kms“^ Mpc“^. 

Previous related works [8, 12] have imposed a hard limit on the width of the potential of 
/ < Mpi, citing both computational reasons (avoiding a divergent direction for the MGMC 
chains to reach convergence) and theoretical ones, snch as inaccnracy of the described po¬ 
tential for / > Mpi and motivations from string theory. Referring to the latter, the PNGB 
potential has been used with the purpose of understanding natural inflation better. One of 
the main points of this approach has focused on the different values of the potential’s width 
//Mpi, in the context of supersymmetry/superstring theories. The overall conclusion favours 
/ < Mpi [12, 20]. However we choose to impose a somewhat weaker prior / < 2Mpi, by which 
time the potential is flat enough that it can commonly generate observables practically in¬ 
distinguishable from those for ACDM. Our motivation to do this is because supersymmetric 
theories have recently become less compelling, given the lack of evidence in their direction, 
and moreover the PNGB potential offers a phenomenological description of dark energy that 
is well behaved, independently of the use of axions during the inflationary epoch. We will 
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Figure 2: Density plots for the equation of state w versus scale factor a of the PNGB model. 
The figure on the left depicts random prior choices from Table 1. On the right, a sample of 
the posterior distribution models obtained from the combined JLA + BAO + HST + Planck 
2015 datasets is shown. The best-fit model is drawn in red, and a zoom is included to show 
the detail of the posterior trajectories given their closeness on the original axis range. 


study the outcome of choosing alternative priors for / later. 

To obtain the observational bounds we utilized the JLA compilation of supernova dis¬ 
tances [21], cosmic microwave background temperature anisotropies, lensing and polarization 
data from the Planck 2015 data release [3, 18], direct constraints on the Hubble constant 
from Hubble Spacee Telescope observations [22], and baryon acoustic oscillations data from 
SDSS [23]. 

4 Results 

The free parameters of our analysis, aside from those standard in any cosmological model, 
are the width / and the amplitude of the potential (2.5). The CAMB code uses dimen¬ 
sionless versions of the field and width parameter, but we will refer to them in their native 
units of mass. The code’s background evolution computation handles all energy-momentum 
components with a normalization equivalent to 87rG/Ocrit/3c^ where pcrit is the present value 
of the critical density. The matter energy constituents are then multiplied using the standard 
density parameter definition Dj = pi/pcrit where the index i is matter or quintessence. This 
means that and are, regarding their units, equivalent to the density parameter for 
the quintessence field 

Figure 2 depicts an ensemble of trajectories for the equation of state w versus scale 
factor a of the PNGB model, drawn from both the prior and posterior distributions. We see 
a very strong tightening of the posterior distribution with respect to the prior, indicating 
that the data are significantly constraining the considered set of models. 

In Fig. 3, a triangular plot with the 68% and 95% confidence contours between the 
amplitude M^, the width //Mpi, the field density parameter the present value of the 
field equation of state wq and the initial displacement of the field cjii/f is displayed. The 
individual posterior distributions of each parameter are also shown. 

The parameter f/Mpi is unconstrained at the upper edge of its prior, and has 68% 
and 95% confidence lower limits of 1.34 and 0.67 respectively. For the 68% and 95% 
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Figure 3: 2D contours of the combined JLA + BAO + HST + Planck 2015 constraints for 
the PNGB model with potential Eq. (2.5). The individual marginalized posterior probability 
distributions of each parameter are also shown. 


confidence ranges are = 0.581 q 26 ~ These probability distributions 

are expected. As / grows the potential flattens, returning the cosmological constant case, 
which makes it impossible to confine this parameter at its upper value, while is sharply 
cut off at the lowest value able to sustain an allowable present density parameter, while fitting 
the data well at that value. Note however that the quoted lower limits on / depend strongly 
on the assumed upper limit of its prior, as models beyond the adopted prior continue to fit 
the data well. We make a detailed analysis of prior dependence in Sec. 6. 

The density parameter is constrained at 68% and 95% confidence as, = 0.69±0.01 
and = 0.69+|5;[j2 respectively. The 95% upper limit on the present equation of state is 
wq = —0.88, whilst the lower cut-off of wq = —1 reaffirms the cosmological constant limit 
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of this theory. This range is as expected for thawing quintessence; while we have permitted 
values of the present equation of state bigger than —1, allowing non-trivial past dynamics for 
dark energy, the cosmological constant case gives a very good fit to the astronomical data. 

Examination of </>i// gives a criterion for the slope of the potential at the start of the 
cosmological evolution. For (j)\/f ~ 0 the feature at the bottom left corner of the (//Mpi, 
(pi/f) duplet shows that only narrow potentials (those corresponding to small /) are allowed. 
In such a steep regime, a slope slightly different from zero would cause the quintessence 
field to evolve too quickly, therefore not reproducing a thawing behaviour. Once the slope 
increases to i^i// ~ 0.3 most of the / range is enabled. Larger values for (j)i/f are favoured 
when / is larger. We hnd a 95% upper limit of (pi/f = 1.9, while there is no lower limit 
since (pi/f = 0 reproduces ACDM precisely. The upper boundary value shows there is a fair 
range of models permitted with an evolving scalar held. However this regime is typically not 
well recovered except for narrow potentials with / ~ 0, because in order to attain such it is 
also necessary for to be very close to Hv, meaning only a narrow sliver of prior space is 
available. The data clearly favours the direction of increasing /. 


5 Thawing quintessence: equation of state parameterization 


An alternative way of parametrizing thawing quintessence behaviour is by use of an approx¬ 
imate analytical expression for its equation of state. Such a study has been carried out in 
this manner in Ref. [15], our main new contribution being the use of more recent data (par¬ 
ticularly Planck 2015) and a comparison with the exact evolution of the PNGB potential. 
The derivation of this solution can be found in Refs. [6, 13, 14], and leads to the equation of 
state 


w{a) 


where 


-l + (l + u;o)a3(-^-i) 


{K - F{a)){F{a) + 1)^ + (R + F{a)){F{a) - 1)^ 

{K - + 1)^ + {K + - 1)^J 

(5.1) 


K 

F{a) 


V 3R(0i) ’ 

“ l)a“3. 


(5.2) 

(5.3) 


The equation of state (5.1) is expressed in terms of three parameters: at present, wq, 

and K which measures the curvature of the scalar held potential at its maximum [6]. For 
K larger than 10 the movement of the held at the start of the evolution is required to be 
very small to avoid a quick roll down. If the held touches the minimum of the potential and 
starts oscillating at a scale factor value near today’s, numerical simulations establish that 
Eq. (5.1) is not valid anymore. In addition to this inaccuracy, for an oscillating potential the 
equation of state would become positive, therefore violating the w < —1/3 condition for a 
dark energy description. For K smaller than 0.5 the held mass becomes very large, implying 
that the Taylor expansion around (p = pi becomes inaccurate because of the rapid variation 
of the held. We have a particular focus on hnding a conhdence range for the curvature of 
the potential K that is better constrained than in past analyses by making use of the exact 
results in Sec. 4. 
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Figure 4: Equation of state versus scale factor for the quintessence field equation of state 
(solid lines) and Eq. (5.1) (dashed lines) with = 0.68. The parameter values in each case 
are (a) //Mpi = 2, cj)-,/f = 1.65, K = 0.98, wq = -0.95; (b) f/Mpi = 1.1, cpi/f = 1.52, K = 
1.02, wo = -0.84; (c) f/Mpi = 0.9, f = 1.43, K = 1.09, wq = -0.77; (d) f/Mpi = 0.5, 
(pi/f = 0.94, K = 1.72, Wo = -0.34. 


In Fig. 4 we compare the behaviour of the thawing equation of state against the nu¬ 
merical solution provided by the quintessence module of CAMB. The approximation works 
very well for tco ~ — 1 and A ~ 1, but becomes less accurate for larger values of wq which 
in turn correspond to smaller values of K. If instead we had adjusted curves generated by 
Eq. (5.1) to the best-fit values of K and wq we would have possibly found a more varied set 
of parameters able to duplicate a broader set of solutions. 

In Fig. 5 the confidence contours for the curvature parameter K are shown. The 95% 
constraint obtained is A = 1.1 ±0.4, where A was calculated as a derived parameter from the 
results of Sec. 4. Given the condition within the CAMB module Quint to attain convergence 
of before allowing the evolution of (pi to commence, the width of the potential dictates the 
value of the curvature. As we discuss in Sec. 6.3, the distribution of A values that emerges 
from simulating the PNGB potential is far from uniform, and hence our constraints appear 
very different from those that sample uniformly in wo and A such as Ref. [15]. In Ref. [15], 
the authors make a distinction between values of A < 1, where they obtain a 95% estimate of 
— 1.21 < Wo < —1.00, and A < 2 which yields —1.14 < wo < —0.92. These values fall largely 
into the phantom region w < —1, and are only partially compatible with ours. However, in 
their case, A is allowed to grow up to a value of 9.95, giving a —2.06 < wq < —1.01. From 
here we can conclude that constraining the cosmological trajectories with a specihc potential 
makes a big difference in preventing this type of models from suffering instabilities. 
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Figure 5: The 68% and 95% confidence level regions for {wo,K) given by Eq. (5.2). 


6 Discussion 

In this section we place the results of Secs. 4 and 5 in context by comparing them to results 
under different choices of datasets and parameter priors, and then compare to constraints 
obtained by previous authors on thawing quintessence. 

6.1 Constraints using Planck 2015 only 

The results of the full-mission Planck observations of temperature and polarization aniso¬ 
tropies of the CMB radiation are widely considered to be the most reliable dataset for con¬ 
straining cosmological models. Moreover, for the standard six-parameter cosmological model 
they are sufficient on their own to fix all the parameters accurately. Here we test whether 
this latter statement remains true with the more general dark energy model by comparing 
results using Planck 2015 only with those of the full dataset described earlier. In Fig. 6 a 
triangular plot with the 68% and 95% confidence contours using the Planck dataset is shown. 
The parameters displayed are the same as those of Fig. 3, but note that in some cases the 
axis ranges have had to be extended. Typically the constraints obtained by this particular 
analysis are considerably weaker than those of Sec. 4. The most striking example of this is 
the width of the potential //Mpp the Planck data alone is unable to constrain this within 
the prior at a confidence range of 95%. A somewhat better constraint is attained for the 
width of the potential M^, where at 68%, = l.l^Qy. This result is linked to the field 

density parameter limit of = 0.64 ^qq^, because of these two quantities’ relation in the 
convergence of the code as mentioned in Sec. 3. There is a noticeable drop in the mean value 
of against the full dataset result of = 0.69; as models with a lower potential amplitude 
are accepted and the field is allowed to roll for a wider variety of initial conditions. 

Most notably, the present value of the equation of state is not well determined, with a 
95% confidence upper limit of wq < 0.27. This result is not very useful for a quintessence 
scenario, given that it does not even require the condition of u^o < —1/3 for a Universe in 
accelerated expansion at present. There is also an unexpected portion of model space around 
wq = -|-1, which is not compatible with broader datasets. We conclude that a CMB-only 
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Figure 6: The 2D contours of the Planck-oT\\y constraints for the PNGB model with potential 
Eq. (2.5). 


analysis of the PNGB model is insufficient to distinguish between the permitted parameter 
combinations, and hence it is necessary to include the geometric data compiled at much 
lower redshift values. This was found also in Ref. [4] for the {wo,Wa) parametrization of dark 
energy. 

6.2 Choice of inverse prior for the PNGB potential width 

In light of the data’s inability to constrain high values of /, we assess the prior dependence 
of our results. We first carried out an analysis which simply extended the upper limit of 
the uniform prior on / to 4 instead of 2, which just has the effect of admitting extra models 
at the higher values of / which make predictions very similar to the cosmological constant. 
As the prior contains more models of this type, constraints on models away from this limit 
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Figure 7: The 2D contours of the combined JLA + BAO + HST + Planck 2015 constraints 
for the PNGB model with potential Eq. (2.5) and parametrization uniform in Mp\/f instead 
of //Mpi. The rest of the parameter definitions are the same as in Fig. 3. 


tighten somewhat, which already indicates a significant prior dependence to any constraints 
which are quoted. 

More importantly, however, there is no particular motivation for choosing a prior uni¬ 
form in /, which did not allow for a clear discrimination between models with a narrow 
width of //Mpi < 1. In order to explore this we changed the parameter with uniform prior 
from f/Mpi to l/(//Mpi). We take the prior on this to be (10“^, 5), leaving the rest of 
our parameters with the same meanings as in Sec. 4. This then means that the cosmological 
constant limit occupies a hnite region of parameter space near 1// ~ 0, rather than the 
potentially-inhnite region when placing a uniform prior on /. Figure 7 shows the results of 
this new choice of prior, where l/(//Mpi) = 2.3 is the upper limit found at a 95% confidence 
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level. This corresponds to a limit on / which is not very different to the one we found in 
Sec. 4, though of course the changed prior modifies the overall shape of the posterior. The 
lower prior limit of 10“^ is close to zero, and therefore indistinguishable by the data due to 
its closeness to the cosmological constant case. 

The rest of the probability distributions are in good agreement as those of Sec. 4, 
indicating that they do not have much sensitivity to the choice of prior. The present equation 
of state wq is a good example of this insensitivity, its limit changing only slightly. However, 
the potential amplitude does become less constrained at its upper values. This is because 
our alternative prior places less prior weight on models which are very cosmological-constant- 
like, hence permitting a somewhat wider range of than our original prior choice. 


6.3 Comparison with previous results 

In this section we compare our results to those of previous authors. Similar studies to Fig. 3 
can be found in Refs. [8, 10-12] and comparison is straightforward after noting a recurrent 
difference in the definition of the variable between those works and ours. As a typical 
example. Fig. 2 in Ref. [8] plots 


1018 GeV’ y^/i/0.65 x 10-3 eVj ’ 

where is the amplitude of the potential, h is the Hubble parameter and / is the width 
of the potential, same as in our case. This choice arises from the fact that the value of the 
critical density is pcrit = (3\//i 10-3eV)"i. The fixed value oih = 0.65 specified by the authors 
gives pcrit = (2.5 X 10-3eV)'i. To translate Eq. (6.1) into our own definition of requires 

= PcritM^ (6.2) 


where fx is the amplitude of the potential in previous references. This choice makes their 
^/{\/h 10-3 eV) ~ 2, which is the minimum amplitude of the potential that would give the 
observed dark energy density. Therefore, « 0.34 in our Fig. 3 is equivalent to ^/{Vh 

10-3 eV) 2.29 in Fig. 2 of Ref. [8]. This somewhat larger is expected, as Ref. [8] 

imposes a hard limit > 0.7 on its models, whereas our data combination constraints prefer 
a smaller A very similar description can be made about the results of the rest of the 
aforementioned references. 

Our analysis results in a significant improvement on the constraint on in relation to 
Ref. [8]; however, the restriction in the same work to wq < —0.965, citing the inability of the 
data to provide a better constraint on this parameter, is invalidated as a significantly larger 
range for the present value of the equation of state is shown in Fig. 3. Concerning the slope 
of the scalar field, the allowed amount of rolling in is in very good agreement with this 
reference. 

Regarding our results for the equation of state parameterization scheme, a comparison 
can be made with Fig. 7 in Ref. [15]. These authors also explore models with wq < — 1; 
however, because quintessence theories do not extend to the phantom domain, we excluded 
that possibility in the prior of our parameter space from the start. Another noticeable 
difference is their use of two datasets and their respective confidence levels; one using BOSS 
data and one without it. When adding BOSS they obtain a confidence contour that is entirely 
in the region wq < —1, but when it is omitted, a small area with rcg > — 1 is still allowed. 
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The shape of the confidence region in Ref. [15] differs greatly from the result displayed 
in our Fig. 5. This is because that article adopted uniform priors in the equation of state 
parameters wq and K. Here, instead, these parameters are derived from a sampling based 
on assumption of the underlying PNGB potential. This induces a prior on those parameters, 
and particularly on K, which is very far from uniform, i.e. the PNGB model realises a very 
different model ensemble from that assumed in Ref. [15]. A simple consequence of using an 
underlying potential is that rco < — 1 is not permitted, but the induced prior on K also leads 
to a substantially different allowed region. Hence caution is required in using the equation of 
state approach to assess the viability of explicit thawing quintessence models such as PGNB. 

7 Conclusions 

In this paper we obtained a probability distribution for PNGB quintessence model parame¬ 
ters, which corresponds to the thawing type. Our analysis was carried out using the codes 
CosmoMG/GAMB, version July 2015. The constraints on the amplitude M^, initial slope 
of the field 4>i/f, and width f/Mp\ of the potential Eq. (2.5) show good agreement with 
earlier analyses, with some reduction in the permitted range of values as compared to them. 
Our results show the continued viability of the AGDM scheme, while showing the extent to 
which models with a present equation of state value larger than wq = —1 and a field evolv¬ 
ing away from the top of the potential (2.5), corresponding to thawing quintessence, remain 
acceptable. 

We then studied the approximate equation of state (5.1), which has been applied in 
Refs. [13, 15]. Under the PNGB assumption, we find much tighter constraints on the cur¬ 
vature parameter K than are found when uniform priors are adopted on the equation of 
state parameters. This difference highlights a strong ongoing prior dependence from the way 
thawing dark energy is modelled, which current data are not strong enough to override. We 
explored the parameter space with different dataset combinations to test their effectiveness. 
We discussed the advantages of a different prior choice and compared with the standard 
parametrization. 

Overall we conclude that the current data are indeed able to meaningfully constrain the 
PNGB model, and in particular force its behaviour to be close to the cosmological constant 
limit. However, in detail there remains a significant prior dependence on the constraints 
obtained, as highlighted by switching the prior from being uniform on / to being uniform 
on 1//. In absence of any clear theoretical guidance on the appropriate form of prior, this 
dependence needs to be kept in mind in interpreting constraints. Hopefully future data will 
emerge with sufficient strength to overcome this prior dependence. 
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